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Abstract 

We develop a low energy effective field theory of a mixture of two species of pseudospin-i 
atoms with interspecies spin-exchange, in addition to density-density interaction. This approach 
is beyond the single orbital-mode approximation. In a wide parameter regime, it indicates the 
existence of the four elementary excitations, especially a gapped mode due to interspecies spin- 
exchange. On the other hand, the spectrum of the effective spin Hamiltonian yielded by the single 
mode approximation can be obtained by quantizing the homogeneous excitation, which is spin 
excitation and is the long-wavelength limit of the gapped mode of elementary excitations. These 
low energy excitations can be experimentally measured by using Bragg spectroscopy. 
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I. INTRODUCTION 



Elementary excitations or collective modes are key properties of Bose-Einstein condensa- 
tion (BEC), and serve as probes of the ground states. The Bogoliubov theory of elementary 
excitations of BEC gives an elegant description of the Goldstone modes associated with the 
spontaneous breaking of U(l) symmetry [l-3]. In recent years, BEC of ultra-cold dilute 
atomic gases has become one of the most active fields in physics. Among the most inter- 
esting topics are BEC of spinor atomic gases [37], for example, spin-1 and pseudospin-| 
gases |s-12], as well as spinless mixtures 13, 14]. As an extension of this topic, it is interest- 
ing to study spinor mixtures with interspecies spin exchange. It has been theoretically found 
that a mixture of two distinct species of pseudospin-| atoms with interspecies spin-exchange 
interaction exhibits interesting features beyond both spinor gases and a mixture of spinless 
gases, especially, in a broad parameter regime, t he g round state is entangled between the two 
species, rather than BEC of individual species 15l-ll7l|. Also, the approach based on single 
orbital-mode approximation has revealed interesting properties of quantum phase transition 



and many-particle quantum entanglement 



16|4l9| 



We expect our work motivates more investigations along this line of research. Spin- 
exchange scattering between distinguishable atoms has been less studied, perhaps because 
of incomplete information on inter-atomic potential. However, we note that interspecies 
spin-exchange interaction can be significant. There are calculations indicating significant 



spin-exchange scattering lengths between distinguishable atom s [2 0|] . Spin-changing scatter- 
ing between distinguishable atoms has indeed been observed [2l| . Experiments on multi- 
component Bose gases often had atom loss due to spin exchanges jjj. [h]. Significant differ- 
ences between singlet and triplet scattering lengths have been observed in 41 K- 87 Rb, m K- 



87 Rb and 6 Li( 7 Li)- 23 Na mixtures 22], implying significant interspecies spin exchanges. It is 
feasible to experimentally realize the systems studied here. One may use, for example, 85 Rb 
and 87 Rb, or 41 K and 87 Rb, as the two species, and \F = l,m,F = 2) and \F = l,m_F = 1) 



as the two pseudospin states 



In this paper, we treat a mixture of two distinct species of pseudospin-^ Bose gases by 
using a field theory approach beyond single orbital-mode approximation. From the point 
of view of field theory, there are four fields, as there are two species of atoms while each 
atom has two relevant spin states. We shall use the path integral formalism to develop a 



Bogoliubov-like mean field theory, in which each field has a specific value in the ground 
state. Excitations are then calculated as small deviations of the fields from those in the 
ground state. 

Previously, elementary excitations in such a mixture have been studied as fluctuations of 
the single-particle orbital wave functions, and it has been restricted to a special parameter 
point, in which the many-body ground state is the so-called entangled BEC 18j. On the 
other hand, when the atoms are all condensed in the same orbital wave functions, there 
are spin excitations described by the effective spin Hamiltonian [161]. In the approach here, 
the spin excitations are obtained as due to spin flipping of the fields that remain spatially 
homogeneous, while the elementary excitations are plane-wave-like excitations of the phases 
of the fields. The former is the long-wavelength limit of the gapped mode among the elemen- 
tary excitation. Furthermore, the low energy effective theory gives the excitation spectrum 
of the effective spin Hamiltonian that is obtained under single orbital-mode approach. 



II. THE MODEL 

Consider a dilute gas of two species of bosonic atoms, the number of atoms in each species 
is conserved. Each atom possesses an internal degree of freedom represented as a pseudospin 
with ^-component basis states | f) and | l), and can transit between the two. This system 
is described by four interacting boson fields, with the Lagrangian density 

where a = a, b represent the two species and a =t, 4- represents the two basis states of 
pseudospin- ~, ip acT = ^ Qcr (x) and \i aa are the field and the chemical potential corresponding 
to the atoms of species a with pseudospin a, respectively, H is the Hamiltonian density 

t f ^Y7 2 _LT/W, \^ „(°«)|„/. l 2 U/, ,l 2 



a a 



(2) 
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where V = V(x) is the external potential, g / , g , and g e are the interaction strengths 
for intraspecies scattering, interspecies scattering without spin exchange, and interspecies 
spin-exchange scattering respectively, proportional to the corresponding scattering length. 



For pseudospin-1 atom, intraspec.es scattering strengths with and without spin-exchange 
are the same [6]. 

For simplicity, we set V = 0, and also assume g^^ = g^ a ^ = g a , gjff' = g^ = g s and 

(ab) _ (ab) 



9-u' = 9±t' = 9 d sucn ^ na ^ = A*c4 — P-a We can define = (ifo(x), ij;^(x)) T , the 

spin density operator S a i(x) = i J a a (. x ) s aa>' l Pav'(. x ) = ^a s% ^a, (i = x,y,z), where s t = T{/2 
is the single spin operator, Tj being the Pauli matrix. Then the Hamiltonian density can be 
written as 



2m Q ' 2 1 1 2 1 1 2 

-\-2g e ( y S ax Sbx + SaySby) + ^9zS a zSbz, 



(3) 



where ^ = # s + gd, g z = 9s - 9d- 

If c/ e = 0, the system is a mixture without interspecies spin exchange, equivalent to a 
mixture of four scalar Bose gases. Note that intraspecies spin exchange does not change 
the particle number occupying each pseudospin state. The Hamiltonian would possess a 
symmetry of U(l) x (/(l) x U(l) x U(l), corresponding to particle number conservation of 
all the four fields. With g e > 0, the symmetry is lowered to U(l) xU(l)xU (1), corresponding 
to the conservations of N a , N b as well as N a ^ — N a ^ + — = 2S Z [18] . 



III. EFFECTIVE LAGRANGIAN OF LOW ENERGY EXCITATIONS 

We consider the parameter regime g e > g z , with all the other parameters fixed. Other 
parameter regimes are studied elsewhere. As a mean field theory, we suppose that in the 
ground state, each of the four fields has a definite value ip aa = ip^e 1 ^" . Then the spin 
exchange term becomes 2g e ip^.ip^ip b \ip^, cos($°^ — $^ — $^ + 3>&jJ. Minimizing the potential 
part of the Lagrangian requires if)°^ = ipa\ = \J Pa/2, V'&t = V'kj, = \/Pb/2 an d ^at - ^aj, - ^01- + 
®b± = 71 ■ We can arbitrarily choose the phases of the four fields under the above constraint to 
describe a ground state, other choices are equivalent in the sense of spontaneous symmetry 
breaking. Therefore in the ground state, Vat = Vv = i>a = \/pa/2, = = V° = 

yPb/2, where p a = N a /Q is the number density of species a, with fl being the volume of the 
system. The chemical potential is evaluated to be \i a = g a p a + \{.9ab — 9e)P/3, where (3 ^ a. 

We now study the elementary excitations. With a deviation from the mean field value, 
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each field can be written as ip aa = + Cacr)e l * aCT , where ( aa is a small quantity. Therefore 



^ 2 4m, 



Pa 



(V$ aCT ) 2 - p a t atT d& c 



+ y PaP6 COS ($ af - $ a4 - $6 t + $ H ) - 2 E CaaC/V , 

a@,aa 

where higher order and constant terms have been neglected. 

Q 

is a 4 x 4 matrix, with 



(4) 



H a H a b ^ 



H a b H b 



(5) 



Pa 



where a ^ (3, 



9aPa + \9ePp, 9aPa ~ \9ePp 
9apa - 9apa + \VePfi 

( 



(6) 



H a b — Pal 



9s — \g e 9d — \g e 
\9d — \g e g s — \g e 



(7) 



Now consider the vacuum persistence amplitude Z = Yi a ,a f 'D&aa'E > ( aa e l I dt S d3xC , from 
which we obtain an effective Lagrangian as a function of <£> only, after dropping the total 
time derivative of $, which does not affect the equation of motion, and integrating over (, 



9* 



p a p 6 cos($ at - $ a4 - $ 6t + $ H ) 



(8) 



In deriving this formula, we neglect the (VC) 2 and ( 3 , ( 4 terms since only low energy 
dynamics is concerned. This Lagrangian has a cosine term similar that in the sine-Gordon 
model. In 1+1D, this term leads to a solution of topological soliton, which has very nontrivial 
contribution to the phase diagram, as discussed elsewhere. However, in this paper we focus 
on the low energy limit in 3 + ID case, in which the fluctuation of $ at — — $^ + ^ is 
largely suppressed and we can make the approximation cos a; wl - x 2 /2. 

The conjugate relation between the phase $ acr and particle number N aai the conservation 
of N a = Yj<j N aa and the fact that the mass term is proportional to (<& -r — — $^ + &bi) 2 
suggest a transformation 

r = c/$ (9) 



where T = (71, 72, 73, 1a) t , $ = ($af, $ai, $bt, ®bi) T , 



( 



U = 



71 71 

73 73 

1 _i 1 _i 

222 2 



1 _i _i 1 
\ 2 2 22/ 



which is orthogonal, i.e. C/ 1 = U T . Then the effective Lagrangian can be rewritten 

c eff = l -{d t v T )A-\d t T) - ^(vr^M-^vr) - y t gy 



where A 1 = UD T Q 1 DU 1 is symmetric, 



D = D 



Hence 



A = UD- l QD- l U T 



1 ' Pa ^ 

Pa 

p 6 

V p 6 y 



2#a 5U - 9e 
9ab - 9e l9b 














g e rj + + g z g e rj_ 

g e r]_ g e r] + - g z J 



where ri± = U^- ± ^), 



M 



-1 



/ £2, ^ 
0^00 

e+ t- 

\ e- e+ ; 



where f± = \( £b -± £b -), 



^0 ^ 





^0 g e p a pb ) 



From Euler-Lagrange equation 



9t(4^k)+n4^k)- dCeff 



d(d t TT) 



0, 
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we obtain the equation of motion of T, 

d 2 t T - M _1 V 2 r + 2AGT = 0, (17) 



IV. ELEMENTARY EXCITATIONS 

For elementary excitations, as characterized by frequency u> and wave vector k, we seek 
solutions of the form of 

r = r exp[-i(wi-k-r)], (18) 
where T is position independent. Hence we obtain 



{-u 2 + k 2 AM~ l + 2AG)T = 0. 



(19) 



The secular equation det(-u 2 + k 2 AM~ x + 2AG) = gives 



_ 0J 2 I gapgk 2 (9ab-9e)Pbk 2 

m a 2mt 



(g a b-ge)pak 2 

2m a 






_ w 2 _|_ gtPbk 2 
mt, 







1.2 













fc 2 (g e /2+g 2 g-) 



2 + 2g 2 v-papb 



(20) 

where f\ = r] + £ + + / 2 = ?7 + ^_ + 

It is found that the four excitations are given by 



1,11 



Y 

CO 



9aPa | 9bPb ^ // gaPa gbPbx 2 | (ga6 ~ 9e) 2 Papb 

m a m b V rn a m b m a m b 



1 



IIIJV 



Bk 2 + A 2 =f VCA; 4 + Dk 2 + A 4 



where 



A 2 \ 2fPb , Pa-, n I 
= l^e(— + — ) -2ge9z\PaPb, 



Pa Pb 



(21) 

(22) 
(23) 



5 = f(^ + ^), C = ^(^-^)2 + ^2^a£L ^ = Se p aft [j2(A-J!.)(£k_&)_2 M /A + 



£) + 2rf(£ + £)l- 



It can be seen that has a gap A, due to the nonvanishing g e , while the other three 
excitations, as Goldstone modes, are gapless. As k — > 0, 



^"2 (5 -2A^ 



k 



C 
4A 2 " 



(24) 
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1,„ D C 



^-^ + -(5 + — )A; 2 + — ^ (25) 



When p a = p b = p, we have B = \g e {^ a + ^)p, D = -2g e g z (g e - g z )(^- + i)p 3 , A 2 = 
2ge(g e — gz)p 2 ■ Note that all our calculations are under the presumption that g e > g z > 0. 

u 2 may be negative in some cases, which means that the mean-field ground state with 
■0at = ^<4 = \l Pa/^i i>b\ — —^bi = \J Pb/2 is unstable and a phase transition occurs. 
From the secular equation we see that u 2 > is satisfied for any k if and only if the 
matrix A is positive definite, as the matrix M _1 is positive definite while G is semi positive- 
definite. This means g a > 0, g b > 0, 4g a gb > (g a b — ge) 2 and g e > g z . The first three 
conditions can be naturally satisfied. If g e < g z , we have B 2 — C = (g 2 — g 2 ) ^ Q ^ & < and 
2BA 2 — D = 2(<7g — g 2 )(^- + < 0, then u 2 n becomes negative for any k and fluctuations 
will destroy the mean-field ground state to form a new phase. 

The parameter point of g e = g z is a point of quantum phase transition. The gap A 
calculated above vanishes at this point, signalling the inappropriateness of the present mean 
field theory for this phase. Indeed the phase at g e = g z is the so-called entangle BEC 
discussed previously by using the single orbital-mode approximation, in which the two species 
are maximally entangled in their collective spins, and BEC occurs in an interspecies two- 
particle singlet state. At g e = g z , the gap calculated in a single orbital- mode approximation 
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19] . An appropriate field theory for 



does not vanish, but is maximal on the contrary 
this phase is under development. 

We can also obtain the correlation function < T[y^(t, x)j v (0, 0)] >, p,u = 1, 2, 3, 4. What 
interests us most is < T[j^(t, x)j4(0, 0)] >. In momentum space: 

Gi(k,u) = -i < [j A (k,u)-f A (-k,-u)} >= 2 9eV+ 2 9z 

By neglecting the k 4 term in ujjy we obtain 

v ' ' J (2tt) 4 u 2 -uj v + i0+ 

= 9en+ - g z 

4n 2 vr 



d r K {vJ{r 2 - vH 2 )A) (26) 



where r = |x|, v = J y + and K (z) is the Modified Bessel Function of the Second 
Kind which has the following asymptotic behavior: 



K (z) 



— In z z< 1 



l z e~ z z » 1 



From G&(x,t), we can also obtain the correlation function of spin-exchange operator 

(£t(x, 0)5 e (0, 0)) oc ( e -*»(x,o) e t74(o,o)j = e iG 4 (x,o) e -iG 4 (i,o) j (27) 

where 1 is a vector of short-range cut-off length. According to (125]) . With g> e > g> z , 
(«S|(x, 0)<S e (0, 0)) decreases with the increase of r. 



V. HOMOGENEOUS EXCITATIONS 



We now consider a homogeneous excitation T(x, t) = T(t), i.e. fluctuations purely caused 
by spin flipping while the orbital wave functions remain homogeneous. In this case the 
effective Lagrangian is L homo = J d 3 xC homo = n[l(d t r) T A~ 1 (d t r) -g e papblj]- The canonical 
momentum conjugate with 7 M is 

dL 

P»=q- = QA ^», (28) 

(/i, v = 1, 2, 3, 4). The effective Hamiltonian is thus 

H homo = Pvin - L = —p^A^pj, + Qg e PaPbll- (29) 

Quantization of these excitations is done by imposing the commutation relation py^jv] = 
i5 pu . Recalling 7^ = U pp $ p and that the N aa and <3> a(T are conjugated variables with 
i N aa, = i$ a p,a*', we obtain 

Pu = ~U up N p , (30) 
where N p represents Ni = N a ^, N 2 = N a ^, N 3 = and A 4 = N^. That is, 

Pi = -j= 2 (N at + N ai ) = -j= 2 N a , 
P2 = -j= 2 (N bt + N bl ) = -j= 2 N b , 

P3 = -K^t - N a± + ~ N K ) = S z , 

P4 = -|(A at - N ai - A fet + N bi ). (31) 

The effective Hamiltonian (12"9"|) can be solved easily. It can be seen that p±, p 2 and p 3 are 
conserved quantities, because N a , N b and S z are conserved. The effective Hamiltonian can 
be rewritten as 

Tj 1 A (de ~ 9z) 2 , 9eV+ ~ 9z , 2 . n 2 / QrA 

H homo = m E PiA ijPj + mger]+ _ gz) Ps + ~x^P* + ngePaWo (32) 



where p' A = p 4 — - p 3 . pi = —N a /y/2, p 2 = —N b /\/2 and p 3 = —5^ are all conserved, 
while the part depending of p' A and 74 is like the Hamiltonian of a harmonic oscillator. 
Therefore the spectrum of H is 

E homo ( Pl ,p 2 ,p 3 ,n) = ^ Eft^Pi + + ( " + ^ 

= ^0 + -Jr-^^ + {n + l )A (33) 

2n(g e r] + -g z ) 2 



where 



£0 = ^[SaiV a 2 + ^iV b 2 + (<7 a6 - £ e )iV a iV fe ] (34) 



is fixed, A is nothing but the energy gap in ([23"]) . In the ground state, S z = 0, n = 0, thus 
the energy is 

E{S z = 0,n = 0) = E + ~A, (35) 

where A/2 is the zero-point energy of homogeneous fluctuation (that is, the k = part of 
( l2Tj) and ( I2"2"j) ). The excitation energy of the homogeneous excitation is 

E homo {S z , n) - E(S Z = 0, n = 0) = ^f e ~ 9 ^ , S 2 Z + n^2g e (g eV+ - g z )p a p b , (36) 

2il{g e rj+ - g z ) v 

where n = 0, 1, • • •. 

It can be seen that the fourth elementary excitation, discussed in the last Section, reduces 
to the homogeneous excitation as k — > 0. 



VI. SINGLE ORBITAL-MODE APPROXIMATION 

For the ground state of a Bose gas, usually the approximation of the single orbital mode 
works very well and is the common practice. For our system, this approximation means 
that only one orbital mode is contained in each field, that is, ip aa ~ a a (f> acT , a a denotes the 
annihilation operator of the orbital mode function <ft acr . Then the spin operator for species 
a is S a = J <i 3 iS a (x) = a\s aa ia ai and thus the Hamiltonian H = J d 3 x'H(x) becomes, up to 
a constant, 

H s = K e (S ax Sbx 4" S a ySby) + J zSazSbz) (37) 

where K e and J z are effective parameters determined by the interaction strengths as well as 
single particle orbital wave functions and energies. Then the total spin S of the system is 
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conserved. In the the uniform case, V = 0, it can be found that K e = 2g e /Q, J z = 2g z /Q, 
where Q is the volume of the system. 

Under the single orbital- mode approximation, ip aa pa a a <f) aa , the elementary excitation, 
that is, the fluctuated phase factor e** aCT with a wave-like dependence on x and t, can only 
be attributed to the fluctuation of the orbital wave function § aa . This verifies the previ- 
ous treatment of elementary excitation using the Gross-Pitaevskii-like equation governing 
the single-particle orbital wave functions. Only in the long-wavelength limit, the gapped 
elementary excitations reduces the homogeneous excitation. 

For a homogeneous excitation, the phase factor e l * aCT is position independent, and be 
attributed to spin degree of freedom. Hence a homogeneous excitation is a spin excitation, 
with the orbital degree of freedom remaining the same as those in the ground state. As such, 
these excitations should be the same as those of the effective spin Hamiltonian H s obtained 
under single-orbital mode approximation. 

As such, the energy spectrum of a homogeneous excitation ( l36l) can be approximately 
equalized with the spectrum of H s , that is 

E s = E homo = ~ Jz) l, S I + ln^2K e (K eV+ - J z )N a N b (38) 

which, for N a — Nj, — N, reduces to 



E s = ~{K e + J Z )S 2 Z + -npK e (K e - J Z )N. (39) 

We have numerically solved the effective spin Hamiltonian H s and compared the result 
with the the above expression of spectrum fl39|) . As shown in Figs. [1] and [21 they fit very 
well for small n. For a Bose gas in absence of a magnetic field, N is very large while S z is 
very small. Hence in ( 1391) . unless n = 0, the first term is much smaller than the second term. 
Therefore, the low-lying states must be those with a certain small S z and with n — 0. This 
firmly indicates that our field theory and the single orbital-mode approximation fit very well 
for low energy excitations. 

This result also confirms a previous perturbative treatment of the anisotropic coupling 



between the collective spins [12|. Both the unperturbed isotropic Hamiltonian and the 
anisotropic perturbation conserve S z , hence the eigenstate is a superposition of states \S, S z ) 
with a same value of S z and different values of S. The expansion coefficients turn out to 
be the "wave functions" of a harmonic oscillator in coordinate S, also giving the spectrum 

11 



400 




FIG. 1. E as a function of S z with n = 0. x = K e /J z and N is the particle number. The '+'s 
are the numerical solution of H s , and the solid line is the plot of E = (K e + J z )S z /4. They fit 
extremely well. 



39]) . The total spin S is indeed equivalent to 74 = ($ a ^ — $,4 — + $&|)/2, as the angle 



between the collective spins of the two species is just 274 — n 



18|- 



VII. SUMMARY AND DISCUSSION 



We have described the low energy excitations of a mixture of two species of pseudospin- 
| Bose gases with interspecies spin exchanges, which entangles the two species of atoms 
when the system undergoes BEC. From the point of view of quantum field theory, we have 
considered a four-component field with various interactions. We developed a low energy 
effective field theory, which can very well describe various low energy excitations in a unified 
framework. As an interesting generalization of the usual Bogoliubov theory for the present 
multicomponent Bose gas with spin degree of freedom, this theory gives four elementary 
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x 10 4 
15i 

E/Jz 




FIG. 2. E as a function of n with S z = 0. x = iT e / and N is the particle number. The dashed line 
represents the numerical solution of H s and the solid line is the plot of E = ^ny / 2K e (K e — J Z )N . 
Note that the low lying excited states correspond to small values of n, for which the low energy 
field theory and the single orbital mode approximation fit well. 



excitations. The most interesting aspect is the gap in one of the four excitations. On the 
other hand, quantizing homogeneous excitations yields the excitation spectrum which can 
be attributed to spin degree of freedom. Interestingly, this leads to an analytical solution of 
the effective spin Hamiltonian obtained under single orbital-mode approximation. 

Notice that in a realistic system in a trapping potential, there is cut-off of Goldstone 
modes due to the trap, thus the low-energy excitations become discrete collective modes j^, 



23]. 



The elementary excitations or collective modes can be measured by using the Bragg spec- 



troscopy, based on two-photon Bragg scattering 24]. Especially, several modes coexisting 



at a given value of momentum transfer can be excited and measured 

13 



25] . For a trapped 



gas, the collective modes can also be measured by perturbing the trapping potential [26 ]. 
Similarly, the excitations discussed in our work can be experimentally measured by using 
the above method. The homogeneous excitations can also be measured, in a way similar to 



24]. 



the measurement of collective modes in a trap, which are not plane waves 

The gap in a collective mode is a feature nonexisting in the usual mixtures, where the 
particle number of each spin state is conserved 22J. The nonvanishing value of g e , which 
accounts for the gap, as well as g z , which characterizes the difference of scattering lengths 
of like-spin and unlike-spin scattering processes, both originate from the interspecies spin 
exchange interaction. Therefore g e and g z are roughly of the same order of magnitude. Many 
experiments have been carried out to measure this interaction [22] , indicating a considerably 
large value of scattering length, which is about 100a#, where a# is the Bohr radius. Therefore 
we expect that experimentally this system can be realized and that the gapped mode can 
be found. 
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